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Abstract 

We show that the parametric correlations of the conductance peak amph- 
tudes of a chaotic or weakly disordered quantum dot in the Coulomb blockade 
regime become universal upon an appropriate scaling of the parameter. We 
compute the universal forms of this correlator for both cases of conserved and 
broken time reversal symmetry. For a symmetric dot the correlator is inde- 
pendent of the details in each lead such as the number of channels and their 
correlation. We derive a new scaling, which we call the rotation scaling, that 
can be computed directly from the dot's eigenfunction rotation rate or alter- 
natively from the conductance peak heights, and therefore does not require 
knowledge of the spectrum of the dot. The relation of the rotation scaling 
to the level velocity scaling is discussed. The exact analytic form of the con- 
ductance peak correlator is derived at short distances. We also calculate the 
universal distributions of the average level width velocity for various values 
of the scaled parameter. The universality is illustrated in an Anderson model 
of a disordered dot. 

PACS numbers: 73.40.Gk, 05.45. -hb, 73.20.Dx, 72.20.My 
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I. INTRODUCTION 



A quantum dot |I| is an isolated 2-D region of a few /im or less in length to which 
M electrons (typically several hundreds) are confined by an electrostatic potential. The 
dots can be weakly coupled via tunnel barriers to external leads in order to study their 
transport properties. For sufficiently low temperatures the conductance of the dot exhibits 
equally spaced peaks with increasing gate voltage @]. This phenomenon can be explained 
within a model where the leads are assumed to behave like electron reservoirs in thermal 
equilibrium at temperature T and chemical potential Ep, and the dot electrons are treated 
as a free Fermi gas with the Coulomb interaction taken into account only through charging 
effects. In this framework each successive peak corresponds to a single-electron tunneling 
event. For an additional electron to tunnel from the leads into the dot and occupy the 
ground state of the resulting Af + 1-electron gas, it has to overcome the Coulomb repulsion 
energy of the dot electrons. Ep must therefore increase from its value at the last tunneling 
event by ^Ep ^ e^/C + A where C is the capacitance between the dot and the reservoir 
and A = -Eat+i — Ej^ is the spacing between the single-particle energies Ej^. This condition 
provides a set of values {Ep} at which tunneling is possible, resulting in peaks in the 
conductance G{Ep) of the dot. The suppression of tunneling between the peaks by Coulomb 
repulsion is called Coulomb blockade. 

Since the charging energy e^/C ^ A, these peaks are equally spaced as a function of 
Ep (for kT < e^/C) and can be observed by increasing the gate voltage which controls Ep. 
The resulting oscillations in the measured conductance are known as Coulomb blockade 
oscillations. In the regime of isolated resonances <C A where Tx is the decay rate from 
the single-particle level Ex into the leads, and for temperatures kT <C A, the conductance is 
dominated by the contribution of the single level closest to Ep. In the temperature regime 
kT ^ Fa, the width of the peaks is ~ kT independent of the individual decay rates. 

The peak amplitudes are found to fluctuate strongly 0. These fluctuations have been 
explained 0] by a statistical theory based on the assumption that the Hamiltonian of the 
dot can be described by random- matrix theory (RMT) @J^. This assumption is valid 
either in the diffusive regime (for weak disorder) or when the classical dynamics of the 
electron inside the dot is chaotic [0 due to irregularities in the confining potential. The 
full probability distribution of the conductance was derived in this framework and found 
to change dramatically when time reversal symmetry is broken due to the presence of a 
magnetic field [Q. 

Recent experiments have been probing the dot's conductance as a function of the ex- 
ternal magnetic field or a changing shape. Of particular interest are the correlations of 
the conductance peaks at different values of the external parameter (magnetic field, shape, 
etc.). For balhstic cavities (open dots), conductance correlations in a magnetic field were 
studied extensively, both theoretically for a large number of channels using the semiclassical 
approximation §| and the supersymmetry method 0, and experimentally Recently 
these correlations were measured also as a function of the dot's shape W^. Such open dots 



are characterized by many overlapping resonances and exhibit properties that are analo- 
gous to Ericson's fluctuations in nuclear reactions ||12|. However, until recently, very little 
was known about parametric conductance correlations in the Coulomb blockade regime of 
isolated resonances, which is analogous to the neutron resonance regime in the compound 
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nucleus [p^fl. 

In this paper we focus on the parametric correlation function of the conductance peaks 



cg{x - x') = G{x)G{x') , G{x) = ^ (1) 

/^2 _ ^2 

where x is any controllable parameter and G{x) is the conductance peak amplitude for a 
given value x of the parameter. In Ref. [1^ we have proposed and derived for the first time 
the universal form of this correlator in the Coulomb blockade regime by using the framework 



of a Gaussian process (GP) [|T5|,|l6| that generalizes the Gaussian ensembles (GE) to situ- 
ations where the chaotic or disordered Hamiltonian depends on a parameter. A particular 
realization of a GP has been introduced in Ref. to study the statistics of avoided cross- 
ings in chaotic systems. Universal parametric correlation of spectra (e.g. level velocities) 
have been derived in Ref. However, the calculation of the conductance correlations 

requires knowledge of statistical properties of the eigenfunctions. The universality of all 
eigenfunction correlators, which implies in particular the universality of (|I|), was demon- 



strated in [15|. Following our predictions, the conductance correlator has recently been 



measured for the situation of broken time-reversal symmetry where x is a magnetic field ||20 



and was found to be in good agreement with theory \T^\. The conductance distributions for 



one-channel symmetric leads were also measured l2l|j2^ recently and are in accord with the 
theoretical distributions derived in Ref. [Q. 

For a dot with reflection symmetry we derive an even more remarkable result: the con- 
ductance correlator is independent of the details of the channels in each lead such as their 
number and degree of correlations, and it reduces to the resonance width correlator for 
single-channel leads. However, the conductance correlations are sensitive to the absence or 
presence of a magnetic field. To demonstrate the universality of our results we test them in 
an Anderson model of a disordered quantum dot. 

Another problem that we address is the difficulty in performing the usual level velocity 
parameter scaling in experimental situations because of the inaccessibility of the spectrum, 
a Coulomb blockade effect. One can trace the variation of the ground state energy of J\f elec- 
trons in the dot as a function of the parameter by following the corresponding conductance 
peak. However it is not possible to measure the level spacing A which is needed for the level 
velocity scaling. Though A may be estimated theoretically, e.g. from free Fermi gas, its 
value may be modified by the electrons' interaction. We find a new scaling procedure, with 
a scaling factor which depends on the eigenfunctions alone and can be extracted from the 
level width or the conductance peak . This should make it possible to extract the scaling 



directly from the conductance peaks data. The scaling factor is interpreted as the average 
rotation rate of the eigenvectors (with respect to x) or as the rms conductance velocity. 

We remark that the conductance correlator in a symmetric dot is identical to the para- 
metric correlator of the eigenfunction intensity at a fixed spatial point of a chaotic system, 
and therefore can also be measured in microwave cavity experiments [^,^. 

The outline of the paper is as follows. In Section II we briefly review the theory of 
conductance in quantum dots in their Coulomb blockade regime using i?-matrix formalism 
|2^-p6[, or alternatively - resonance theory p6[. Both formalisms are particularly suitable 



for making the connection to RMT In Section HI we derive and demonstrate the 
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universal parametric correlations of the conductance peaks. We consider symmetric dots 
and show that in this case cg{x — x') is not only universal upon an appropriate scaling 
of the external parameter, but also independent of the channel information. We compute 
its universal form using a simple GP for both cases of conserved and broken time reversal 
symmetry. In section IV we introduce the new scaling procedure based on the eigenfunction 
rotation rate. The scaling factor diverges logarithmically for the orthogonal case but we 
describe a way to regularize it. Using perturbation theory, we derive analytically the width 
correlator at short distances and express the rotation scaling in terms of the RMS of the 
width velocity. We also derive a semiclassical expression which is system-dependent for the 
rotation scaling and compute it for the case when the parameter is an external magnetic 
field. This gives us the magnitude of the correlation field. In Section V we compute the 
universal distributions of the average width velocity for various values of the scaled Ax. 
In Section VI we discuss the universal parametric correlations of the conductance peaks in 
the general case of asymmetric dots with symmetric or asymmetric leads, and express the 
rotation scaling in terms of the RMS of the conductance velocity. Finally, our analytical 
results are verified in the completely solvable N = 2 GP (Section VII). 



II. CONDUCTANCE IN QUANTUM DOTS 

A. i?-Matrix Formalism 

It is convenient to discuss the conductance of the dot in the framework of i?-matrix the- 
ory, originally introduced for nuclear reactions that proceed through the compound nucleus 
formation and decay [^,^,0] and more recently applied to quantum dots P, p7| , p8| . We 



briefly review this formalism in the context of quantum dots but for the general case where 
there are several open channels at each lead. Let us consider a planar cavity-like dot in the 
x-y plane with left and right leads of width Di,Dr attached to it at the lines of contact 
X = xi,Xr', < y < Di j. which we denote by C. The dot Hamiltonian has energies Ex 
whose corresponding eigenf unctions '^\{r) vanish at the walls and satisfy an homogeneous 
boundary condition 

d-^x/dn - hi^r'^x = for f e C (2) 

where h is the normal to C (i.e. n = ±x) and hi^r are constants. 

We denote by 0* a complete set of transverse wavefunctions, where c is a channel in- 
dex and i = l,r. For an open channel 4>c{y) = \J2/ Di sin(K*?/), where = cix/Di (c = 
l,2,...,Aj) is the quantized transverse momentum. Inside each lead we can expand the 
scattering solution $ at energy E as 

m = j:<{^my)- (3) 

c=l 

The i?- matrix, which relates the values of du\t/dn — hiiu\t to those of at 

be expressed in terms of the resonances' eigenfunctions and energies through Ri>c'-ic{E) = 

Eylxy'c'x/iEx-E) where 
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ycx = \h- dl^l*^, (4) 



is the reduced partial width amphtude for the decay from level A into channel c in lead i. 

If the resonances are isolated, i.e. their total width (see below) is much smaller than 
their average spacing A, only one resonance whose energy Ex is closest to E contributes to 
the i?-matrix, and the cross section for scattering from channel c in lead i into channel c' in 
lead i' {i 7^ «') is given by the Breit-Wigner resonance formula 



pi pi 

Here r*;^ and Ycx respectively the partial width and partial width amplitude of the 
resonance level A to decay into channel c in lead i 



Kx =1 7i r ; iix = pmvcx (6) 

where k], is the longitudinal channel momentum ^^/c*^/2m + Ti^k^^ /2m = E, and is the 
penetration factor to tunnel through the barrier in channel c. r\ = J2c ^Ix total decay 

width of level A into lead i. 



The conductance of the dot at energy E is given by Landauer's formula 9{E) = 
(e^/h) J2cc' I Src';ic At zcro temperature (or for kT ^ Fa) the conductance as a function 
of the Fermi energy is given by g{E = Ep), while at finite temperature the conductance 
g{E) due to one resonance level A has to be convoluted with the derivative of the Fermi 

function /(e) = (l +6^/^=^)"^ 

G{Ef) = -j dEg{E)f{E - Ep) . (7) 

In the temperature regime F^ ^ kT ^ A where Coulomb-blockade oscillations have been 
observed the peak of g{E) at E = Ex is much narrower than the peak of —f'{E — Ep) 
there and we can approximate (0) by 



G{Ep) = cosh-^ " , " . (8) 

The conductance has therefore a resonance shape as a function of Ep centered at the reso- 
nance energy Ep = Ex with peak amplitude of 

^ e2 TT F'F^. 

G = (9) 

h 2kT T{ + T\ ^ ^ 

and width of kT. This formula has been derived in using a different method. 

We remark that a different modelling of a quantum dot assumes point-like contacts where 
each lead is composed of several such point contacts. One arrives at the same expression 



for the conductance [30|, except that each point contact is considered as one channel 



and the corresponding reduced partial width is |]3l[ 7ca = JcucAA/ir'^xi^c), where A is the 



area of the dot, A is the mean level spacing and etc is a couphng parameter of the dot to 
the lead that represents a tunneling probability through the barrier. 
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B. Resonance Theory 



Eqs. (^) and which relate the partial widths to the resonance eigenfunction in the 
dot, play a major role in the statistical theory for the conductance peaks (|]), and it is 
therefore instructive to derive them directly from resonance theory p6| without using the 
i?-matrix formalism. 

A decaying state is described by the condition that the wavefunction has only an outgoing 
component in the external region (instead of the boundary condition (0) used in the i?-matrix 
formalism). In our case a general solution $ at energy E has inside the leads the form 
with Uc{x) = AcU'^{x) + BcU^{x) where u^{x) are incoming and outgoing longitudinal waves 
propagating through a potential barrier Vc{x) in channel c. The conditions i?c = in all 
channels determine a discrete set of complex energies = E\ + 5E\ — iT\/2, where 5E\ 
is a real shift in the resonance energy and T\ = J^c^cx is the total resonance width. An 
equivalent boundary condition for determining the complex resonance energies is that the 
logarithmic derivative at the interface of the leads with the dot is equal to the logarithmic 
derivative of the outgoing wave 



^1 M 



d 

/i+ = — In I 

~ dn 



(10) 



are generally complex, e.g. hf^ = ±ikl in the absence of barrier (when is a plane 
wave propagating in a direction away from the dot). The modification of the boundary 
conditions from (^) (where hi are real) to (|10|) allows the wavefunction to leak into the 
leads. The partial decay rate Tc\/ti into channel c of lead i is given by the current through 
the cross-section of the lead 



j_cX 

h 



h 



2im 



dl m 



dn 



dn 



h 



2im 



.duc 
dn 



Ur 



dK 
dn 



(11) 



Using the boundary conditions ( p^ we find 



cX 



m 



(12) 



which is exactly Eqs. (13) and (H) provided the penetration factor Pc is defined by 



(13) 



In the absence of a barrier = 1, and we expect P^ ^ 1 in the presence of a barrier. 



C. Statistical Model for Dots with Multi-channel Leads 

Eq. (|^) indicates that the irregular fiuctuations of the conductance peaks arise from 
fiuctuations of the level widths F*^. We assume that the channel penetration factors have 
negligible dependence on the parameter, and in any case they vary smoothly with x. These 
fiuctuations must then come from the eigenf unctions \I'a inside the dot, and have recently 
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been accounted for by a statistical theory based on the assumption that the Hamiltonian 
of the dot can be described by a x random matrix H taken from the appropriate GE 
0. The dot's eigenfunctions \1/a are expanded in a complete set of solutions inside the 

dot with the given energy Ex, ^I^a(^) = V'^P^lr); ^^"^ expansion coefficients ip^ are 

assumed to constitute the components of the A-th eigenvector of the random matrix H. 
The partial width amplitude to decay into channel c can then be expressed as a scalar 
product 

TV 

llx = E = i€ I V^a) , (14) 



where 0* is the channel vector defined by 0*^ = ^c-^c/"^ / d-^ 0c*P/i- ^^e point con- 



tact model, a similar formula (14) applies but with 0c^t = y <^c-^^/^P)l(rc)- Thus in the 
corresponding iV-dimensional space, the partial width amplitude of a level is simply the 
projection of its corresponding eigenstate ipx on the channel vector 0c- We remark that the 
scalar product in Eq. (0) (that will be used throughout this paper) is different from the 
original scalar product defined in the spatial region occupied by the dot. 

We define the channel correlation matrix M of lead i by M*^, = 7" 7c/- Using the scalar 
product expression ([M]) for the partial width amplitude and the GE average ^f^^^i = -^^iifi', 
we can rewrite M in terms of the scalar product of the corresponding channels 

Ml, = ^(0^ I 0^) - (15) 

The norm of a channel vector is determined by the corresponding mean partial width through 
(0c I 0c) = . Channels c, c' are called uncorrected if 7c*7c' = and equivalent if 
r* = r*,. In general channels can be correlated or non-equivalent and the channel vectors 
0c can therefore be non-orthogonal and have different norms. Channels in different leads 
are assumed to be uncorrelated. The matrices M^'^, thus contain all the channel information 
relevant for the statistical description of the dot's conductance through Eqs. (||), (||) and 
(0)- 



III. CONDUCTANCE CORRELATIONS IN SYMMETRIC DOTS 

In experiments on quantum dots one can vary an external parameter x such as the shape 
of the dot or the strength of an applied magnetic field and trace the corresponding change 
in the amplitude of a given conductance peak. We are then interested in calculating the 
conductance peak correlator (0), assuming the electron's dynamics inside the dot is chaotic 
(or disordered) at all values of x. In this section we show that the conductance correlator is 
universal upon an appropriate scaling of the parameter x and we obtain its universal form. 
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A. Gaussian Random-Matrix Process 



To calculate the correlation between conductance peaks that belong to different param- 
eter values it is necessary to use a random-matrix model which incorporates a parametric 
dependence. Such a model, the Gaussian random-matrix process (GP) which constitutes a 
natural generalization of Dyson's Gaussian ensemble (GE) |^ has recently been proposed 
and applied to the calculation of universal parametric correlation functions in chaotic and 



disordered systems |T5 



A GP is a set of TV X random matrices H{x) whose elements satisfy 



HxJx) = 



Hxa{x)H,,{x') = ^f{x - x')gfl^^ , (16) 

where g^xa~i^u — ^X/i^cru + ^Xu^afi and g^xa~i^u — '^^Xu^a^i- For situations with time-reversal invari- 
ance we have /3 = 1 and H{x) are real symmetric matrices, whereas for broken time-reversal 
symmetry (3 = 2 and H{x) are complex Hermitian. The corresponding GPs are termed 
Gaussian orthogonal process (GOP) and Gaussian unitary process (GUP), respectively. 
A GP is characterized by the short distance behavior of its correlation function / 

/(x - x') ~ 1 - K I X - x' |\ (17) 

where k and (0 < ?7 < 2) are constants. The GP unfolded energies = Ex/ A (where A 
is the average level spacing) satisfy a diffusion law at short distances 



{Asxy = D\Ax\'^+0 (l Ax p") (18) 

where 



^ (Aex^ 4:Nk 

D = hm ^ = -— 19 

plays the role of the diffusion constant. Upon the parameter scaling 

x^x = D^'^x (20) 

the correlation function /(x — x') (^) becomes / ~ 1 — | x — x' |^, and all correlators 
become universal functions of | x — x' |''. This universality can also be demonstrated through 
Dyson's Brownian motion model P3|. The t] = 2 GPs are the only class of processes that 
are continuously different iable in the parameter and therefore the ones suitable to describe 
most physical systems [Q. For rj = 2 the diffusion scaling (|20|) reduces to the level velocity 

scaling D = C(0) = (dex/dx)'^, first introduced in |]18 . 
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B. The Conductance Parametric Correlator 



When the dot is symmetric under reflection around the y-axis we have \1/a(— x,?/) = 
±\E'a(2^5 y), thus r^_;^ = T^^^ = Tcx and the peak conductance (H) becomes 

We take the dot Hamiltonian H{x) to be a member of the appropriate GP with energies 
E\[x) and eigenf unctions ip\{x). The conductance correlator (|T]) reduces in this case to the 
level width correlator 

cg{x - x') = cr(x - x') = f{x)f{x>) , f (x) = ^^S^ (22) 



A 

where r(x) = Yl, I (0c | '^'(x)) P (see ([l^ ). Here and in the following we omit the eigenvector 

c=l 

label A and the lead label i. It is clear from the discussion in Section III. A that this correlator 



becomes universal under the diffusion scaling (|20|), for any given set of channel vectors (pc- 
Furthermore, due to the invariance of the GOP (GUP) under an x-independent orthogonal 
(unitary) transformation, this correlator can depend only on the correlation matrix M in 
(p!5|). However, we will prove below the stronger result that the width correlator is also 
independent of M. 

Since M is Hermitian and positive-definite, it can be diagonalized by a unitary trans- 
formation U under which the channel vectors (pc transform into an orthogonal set of eigen- 
channels (pc {(pc = J2c' ^c'Uc'c) and 

Mcc' = ^{^c\ 0c') = wXc' (23) 

where w'^ = (0c I 4>c)/N are the (positive) eigenvalues of M. Since the total width (in each 
lead) is invariant under a unitary transformation, we have F = J2c I 7c P= J2c I 7c P (where 
7c = (0c I V'a) are the partial widths to decay to the eigenchannels) and the width correlator 
thus depends only on the eigenvalues if^ of the correlation matrix. Defining the normalized 

eigenchannels | 0c) = (v^Wc) | 0c) we can express the width in terms of the partial 
widths to decay into an orthonormal set of channels 

A 

F(x) = Nj2wlt,{x) , fc(x) =1 (0c I ^(x)) I' , (24) 

c=l 

and the numerator of the width correlator (123) becomes 



F(x)F(xO - f' = J2wt ( fc(x)fc(x') - fc' ) + E ^cV' ( fc(x)f c'(xO - f c fc' ) . (25) 

c^c' 



We now observe that the averaged quantities in ( p5D do not depend on the subscripts c, c' 



but only on the orthogonality relation (0c | 0c') = 0. Indeed, if we chose a different pair 
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(pci <|3c' satisfying ((^c | <^c') = there would exist a unitary transformation U such that 



U 



c) and I (pc') 



U 



>c') which could be used to rotate H{x) at each x 
into UH{x)W with eigenf unctions U \ ip\{x)), leaving Tc{x), Tc'i^') unchanged. Since the 
GP probability measure P [H{x)] D[H{x)] (see [|l6l) is invariant under global rotation, the 
corresponding averages remain unchanged as well. We then rewrite (E^) as 



T{x)T{x') 



E 



T,{x)T^{x')-T^ 



+ 1 E 



ri(x)r2(x')-ri 



, (26) 



where the subscripts 1, 2 refer to any orthonormal pair 0c, <Pc' (c 7^ c'). Now the cross-channel 
correlator (second correlator on the r.h.s. of (|26|)) is related to the autochannel correlator 
(first term there) by 



ri(a;)r2(x')-ri 



N -I 



ri(x)ri(x')-ri 



To prove ( P7| ) we use a complete basis of orthonormal channels 



N 



plying the two normalization conditions Yl, I 7i(^) 

i=l 

we get 



N 

P= 1 and E I 1 2= 



(27) 

1,...,A^. Muhi- 
1 and averaging, 



A^l 71 (x) |2| 71 (xO |2 + N{N - 1)1 7i(x) |2| ^^{x') 



- 1 . (28) 

Relation (p8[) leads immediately to (pT]) . According to (pT]) the cross-channel correlator is 
negligible in the limit N 00 compared with the autochannel correlator. From ( p6|) we 
obtain 



T{x)T{x') 
and therefore 



E 



N 



X 



1 1 
Wc^d 



ri(a;)ri(x')-ri 



ri(x)ri(a;') . 



(29) 



(30) 



This result ( pO]) is remarkable in that the conductance correlator in a symmetric dot 
is not only universal but also independent of the number of channels in the leads, their 
associated mean decay widths and the correlations among them; it is given by the correlator 
of a dot with single-channel leads. The universal form of this function can be computed 
using any GP [IT3. We choose the simple process pi 



H(x) = Hi cos X + H2 sin x 



(31) 



which is an ?7 = 2 GP with a correlation function f{x — x') = cos{x — x'). Hi,H2 are 
independent matrices belonging to the appropriate GE. Fig. shows the results for both 
the GOP and GUP. We find that cr{x — x') is approximated very well by 



crix 



x') 



1 + (x — x'Y/a'^ 



(32) 
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namely a Lorentzian in the orthogonal case and a squared Lorentzian in the unitary case. 
Using 300 simulations with = 150 and estimating the statistical error at each x, a least- 
squares fit to (^) gives «! = 0.48 ± 0.04 and a2 = 0.64 ± 0.04. The numerical simulations 
of the correlators (diamonds) and their best fit to (0) (dashed) are shown in Fig. |l|. We 
have also confirmed in our simulations the independence of the details of the channels by 
computing the correlator (p2D with various sets of channel vectors (pc and verifying the 
agreement with (P^. Comparing the orthogonal and unitary results shows that breaking 



time-reversal symmetry accelerates the decorrelation. In Section IV. C we use perturbation 
theory to derive an exact analytic expression for cr to leading order in Ax. This leading 
order differs from that of (|32|), and thus expression ( |32D is just an approximation. 

To test the universality of the conductance correlator, we model a disordered dot by the 
two-dimensional Anderson Hamiltonian [ p5| , p^ 

H = ^(wa + w„)aj,a„ - (33) 

a (a/3> 

where the second sum is over nearest-neighbors lattice sites only. The on-site energies Wa 
are distributed uniformly over the interval [—W/2, W/2]. The Ua describes the non-random 
component of the potential whereas the phases 6ai3 are determined by the magnetic flux. 
We use a Ux y. Uy = 27 X 27 lattice folded into a cylinder whose symmetry axis lies along 
X, imposing appropriate boundary conditions. The left and right leads are represented by 
arrays A'-'^ of mi x m2 lattice points with total widths given hy T\ = J2 ("^q)^ I ^A('^a) P 

such that each point contact constitutes a channel with coupling (f^)^ [pO| . 

We introduce a parametric dependence by applying a step potential along x, u{x,y) = 
UqQ{x), using its strength uq as the parameter. The results for the conductance correlator, 
which are presented in Fig. |I| (left) for 2x2 leads and = 1, are in excellent agreement 
with the GOP prediction (with p = 1). 

We break time-reversal symmetry by applying a magnetic field B pointing along the 
cylinder axis, corresponding to 9ai3 = ^$/$o in (^) if ^a,^/? are neighbors along y and 
0af3 = otherwise, where $/$o is the magnetic flux associated with B in units of the flux 
quantum $o = hc/e. We take $/$o = 1/4 to achieve maximal symmetry breaking. The 
parametric dependence is introduced through a step potential and the results are presented 
in Fig. ^ (right) using the same lead arrangements as above. The agreement with the GUP 
prediction ( P^ (with /3 = 2) is again very good. 

We remark that the partial width is analogous to the eigenfunction intensity at a fixed 
spatial point f. It follows that the parametric eigenfunction intensity correlator in a chaotic 
system is identical to the width correlator cr (0)- Denoting by \l/^(r^ the eigenfunction at 
r for a value x of an external parameter, we have 



I{x)I{x') = crix - x') , J(x) =1 ^^(f) r • (34) 



The eigenfunction intensities are measurable in microwave cavity experiments [22,23| and 
our predictions for cr{x — x') can be tested there as well. 
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C. Relation to the Parametric Overlap Correlator 



It is interesting to note that the result (|32|) is identical (within statistical errors) to our 
result for the universal form of the averaged parametric overlap [ffSl 



oix — X 



{^x{x) I Va(x')) P 



(35) 



We now derive an exact relation between the width and overlap correlators (for finite A^) 
and show that they are identical in the limit oo. Choosing a complete set of vectors 

(pi ] i = 1, . . . , N the overlap correlator can be written as 



o[x 



X 



n I 



i+3 



(36) 



Since the various correlators in (|36D are independent of the particular channel i or the pair 
z, j of orthogonal channels we can rewrite it as 



o\x 



X 



iv ri(x)ri(x')-ri 



A^(A^ 



7i*(x)72(x)7i(a;')72*(a;') + c.c.J + - . (37) 



To proceed we need to use a relation between the two parametric correlators on the r.h.s. 
of (P7D, the width correlator and the "exchange" correlator: 



7i*(a;)72(x)7i(a;')72(a;') + cc. 



15 N 
2N-1 



r,(x)r,(x')-r. 



-21 



(38) 



At X = x' , this relation is easily derived using the results of Appendix A, and for small 
X — x' 7^ 0, it can be proven perturbatively (see Appendix B). At the end of this section we 
provide a proof of this non-trivial relation for any x, x' . Using ( p8D in ( |37D we find 



o{x 



X 



N 



pN + 2 



ri(x)ri(x')-ri 



+ 



N 



(39) 



Finally, using the GE relation tj-ti = 2{N - l)/N'^(pN + 2) (see Appendix A), we find 
the exact relation between the overlap and width correlator 



o[x 



X 



1 



1 77 I cr(x — x') + — . 



N 



(40) 



In the limit N ^ oo these two correlators become identical. 

We now return to prove ( pSD . We first rederive the independence of the width correlator 
of the matrix M in a different way. The channel vectors 0c can be transformed by a linear 



transformation 



J2c' (f'c'Fc'c to orthonormal channels 



') = Sec'- In terms of F, the 



original channel matrix is {(pc \ 4>c') = {F'^F)cc'- A simple solution for F is F = {NWy/^U 
where U is the unitary matrix that diagonalizes M and W is the diagonal matrix composed 
of the eigenvalues w"^ of M (see previous section). In this case FF"^ = NW is diagonal. 
However, we can take advantage of the freedom to make a unitary transformation R on 
the orthonormal eigenchannels of the previous section to obtain another orthonormal set 
0c- This leads to F = R{NWy^'^U which in general has a non-diagonal FF^ = NRWR^. 
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Calculating the width correlator directly from the expression r = ^^FF^ for the total 
width, we now find 



r(x)r(x') - = ^ Uff^), 



ri(a;)ri(a;')-ri 



+ \j:iFF^UFF^). 



7i*(x)72(x)7i(x')72(a;') + c.c. 



(41) 



The third term on the r.h.s. of (|41|)) is just the cross-channel correlator and according to 
(^) can be combined with the autochannel width correlator (the first term on the r.h.s. of 
(PI)). Notice that for the choice made of F in the previous section (FF^ = W diagonal) 
the second term on the r.h.s. of (^T]) vanishes, unlike the present case. However, since one 
must recover the same width correlator in ( |2^ ) and in (^), we conclude that the first two 
correlators on the r.h.s. of (^If) must be related. Indeed, if we assume relation ( pS]) to hold, 
the two terms in (^Tj) combine to give 



T(x)T(x') - r 



Tr(FF^)2 - 



A^- 1 



(TrFF"f)2 - Ti{FF^y 



ri(a;)ri(a;')-ri 



• (42) 



Since Tr(FFt)2 = TrM^ = Ec^^t and Tr(FFt) = TrM , Eq. (j^ is then identical to 



IV. THE ROTATION SCALING 

Although all parametric correlators become universal upon the diffusion scaling (pOD, per- 
forming this scaling is not always feasible. In particular, the scaling factor D is not directly 
accessible in quantum dot experiments. As a result of the Coulomb blockade suppression 
of tunneling, the energy levels probed as the gate voltage is varied are the ground states 
for successive electron numbers rather than the excited states for a fixed electron number. 
Therefore one can measure the charging energy e^/C but not the level spacing A which is 
necessary to determine D. While in principle it is possible to measure the variation of the 
ground state energy of the dot as a function of x by tracing the peak amplitude, the mean 
level spacing A can only be estimated theoretically from e.g. a Fermi gas model. Since the 
interaction and the dot's shape affect A, it is difficult to extract accurate values for D from 
the data. 

In this section we introduce a new scaling which is calculated from the eigenfunctions 
alone, or alternatively from the level widths. This makes it possible to extract the scaled 
parameter directly from the conductance data (see also Sec. VI. C), so that our predictions 
for cg{x — x') could be tested experimentally. 



A. Eigenfunction Rotation Rate and Parameter Scaling 

We consider the quantity | A'jx p = | 'Jx{x') — ■yx{x) p where 7a (a;) = (0 | "ipxi^)) is the 
partial level width and is a channel vector which we take for convenience to be normalized. 
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This quantity is independent of our choice for since the GOP (GUP) is invariant under 
an x-independent orthogonal (unitary) transformation that will transform into another 
normalized vector. From first-order perturbation theory we have 



A7A 



E 



{iP^jx) I H{x') I i,x{x)) 
Ex(x)-E„(x) 



7^(x) +iuJxix)-fxix)Ax , 



(43) 



where ujx{x) is an arbitrary real function that corresponds to the freedom we have in choosing 
the phase of ipx. We first make the usual choice ux = 0, but we shall see below that another 
choice would be useful for computational purposes. The variance of A7A is calculated in 
two steps. First, we average over H{x') keeping H{x) fixed using the conditional two-matrix 



distribution 15 



P [H{x') I H{x)] = P [H{x), H{x')] I P \E{x)\ 



oc exp 



/5 



{Yx\H{x')~]H{x)\ 



with f = f (x — x') . We can write the conditional averages as 



fH,x)iK,-fH^,) = ^9ltvi'^ 



(44) 



/iA 



(45) 
H' I 



Here and in the following we use the notation Hf^x = {i^tj, \ H \ ipx) and H'^^^ 
i/jx), where matrix elements are calculated in the basis of the eigenstates of H{x). Primed 
quantities take their value at x' whereas non-primed quantities take their value at x. To do 
the conditional averaging we write | A7A P in the form 



r= E 



fH^x) {H'*x - fH^ 



'uX) 



{Ex-E,){Ex-E, 



which is convenient since P [H' \ H] depends on H' — fH (note that H^x 
averaging over H' we get 



A7A 



^(i-/')E 



7m 



2/3 



(46) 
0). After 

(47) 



Next we perform the average over H. Since the joint eigenvalue and eigenfunction distri- 
bution factorizes in the GE and | 7^ P is /i-independent we have, using the short-distance 
expansion of f{x — x') 



where 



A7a |2 = (3\ 7a \mZK I Ax 



1 



(48) 



(49) 
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Dyson |^ obtained an exact result for a related quantity, 

However, i?A in (^91) is near the middle of the spectrum where the level density is roughly 
constant whereas in (|50|) it ranges over the whole semi-circle, hence Z ^ Z^. It is more 
suitable for our purpose to unfold the levels and approximate the sum in ( ^Q]) using the 
two-point cluster function 12(e) 



00 



The integral diverges for /3 = 1, as does Dyson's result (50). For /9 = 2 it can be calculated 

1 

using the simple Fourier representation Y2{e) = sin^(7r£:)/(7r£)^ = / dke~'^'"^'^^ {1— \ k |), 

-1 

giving Z = 2/3 compared to Dyson's Z^ = 1/2. 

Eq. (ESI) can be rewritten (for Ax 0) as an amplitude diffusion law 



\A^x\^/\lx\^=f3R\Ax\^+... (52) 
in analogy with the level diffusion law ([18|) , where 



Ax- 



1 


A7A 


2 




2 


Ax 


V 



R = lim / ; . (53) 



In the GP we have R = kNZ. Eq. (^) motivates the parameter scaling 

X ^ X, = i?^/^x (54) 

which we call the rotation scaling as explained below. Under this scaling, just like the 
diffusion scaling (|20|) , the process correlation function ([T7| ) becomes 

/ ~ 1 ~ \ Xr — \^ , (55) 



we 



and since all correlators depend on N{1 — /) rather than on N and / separately []15 
conclude that all correlators are universal functions of I x,. — x' I''. We remark that R can 



also be written as R = f3 lim | A7A p/ | Ax |'' where 7a = 7a/v I 7a P is an "unfolded" 

Ax^O 

amplitude in analogy with the unfolded level used in the definition (|l^). 

For the case of broken time- reversal symmetry (GUP), the amplitude 7a is generally 
complex when evaluated according to (^3]) with ux = 0. Note, however, that the eigenfunc- 
tions ipx{x) are determined up to an (x-dependent) phase, and we can choose this phase 
such that the amplitude 7a = (0 | "^a) is real for any x. For this choice to be consistent we 
must fix co'a(x) such that a;A(x)7A(x)Ax cancels the imaginary part of the sum on the r.h.s. 
of (PD, giving 
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A7A = E (56) 



which is real as well. Since only the real part of the matrix elements of H is used in (p6D, 
the expression (E3|) for R is now replaced by 



fl.Hn, 1(424!. (57) 

In practice it is more convenient to use this definition, according to which 7^ = y/Tx and is 
thus real. 

In most physical applications, the dependence of if on a; is analytic and the GP must 
be differentiable, i.e r] = 2 [jl6|. The scaling factor is then given by 



R=—\^?, (58) 

which can be interpreted as the eigenfunction rotation rate. To see this we note that as x 
increases, the eigenfunctions rotate in the Hilbert space and this rotation is characterized 
by an anti-Hermitian operator Q{x) 

^\^x) = Q\ V-a) (59) 

whose matrix elements in the basis ip\ are given by fi^A = (V'/i I dipx/dx) = ( 1 — (5^a) ^^^x-e^^ 
(see (|4^)). The basis-independent quantity iV~^Tr(f2tJ7) is obtained from (|59|) with the help 
of a complete set of orthonormal vectors 0^ 



-Trinm) = - E I I d^Px/dx) \^ = R (60) 



Xfj, 



exploiting the (/)^-independence of the average in the sum in (^) and the fact that 7^ = 
R is therefore the variance of the matrix elements of the scaled rotation rate matrix 



'N^l{x) = Vt[x)/^ \ 7a in analogy with D being the variance of the scaled level velocity 
{dEx/dx)/A. It also follows from (|5|) that R is the norm of the eigenfunction derivative. 



R = {diPx/dx I dipx/dx) . (61) 



If the system obeys the GE statistics at each x, then d'yx/dx = and the quantity 



I d^x/dx P in the definition of R gives the variance of d'jx/dx. This should be contrasted 
with the level velocity scaling D, where the physical system can have a drift such that 
dex/dx 7^ which has to be removed. 

While the values of the scaling factors D and R are system-dependent, the ratio R/D is 
universal and is given by 

R/D = ^TT^Z . (62) 
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For the GUP Z = 2/3 and 



R/D = — {GUP) . (63) 
3 

If we use Dyson's Zj:, in (pHD, we obtain R/D = 7r/2. The resuhs from simulations are in 
clear agreement with R/D = 7r^/3, as can be seen in Fig. 

It is important to note that R diverges in the orthogonal case (in the limit Ax 0) 
because of the divergence of Z in (0,^) and the rotation scaling, as defined above, is valid 
for the unitary case only. However, it is possible to fix this flaw as we discuss in the next 
section. 



B. Rotation Scaling in the GOP 

We now analyze the divergence of R in the orthogonal case and show that it is nevertheless 
possible to use it to perform the scaling. We note that | A'jx I 7a I ^ is bounded from above 
by 4 and is therefore finite for any Ax. This is in contradiction to (|52|) , according to which 
this quantity diverges for any fixed Ax > 0. This indicates that while the perturbation 
expansion (|43|) is valid for a given member of the ensemble (for Ax sufficiently small), the 
ensemble average of the leading order term diverges due to contributions from members with 
small level spacings. Nevertheless, since | A'jx P/ I 7a P is finite, we can define a scaling 
factor R{Ax) (in analogy with (^)) for any finite value of Ax 7^ by 



^(Ax) = 



1 1 


A7A 


2 


P\lx\ 


2 


1 Ax 1 





(64) 



R diverges (in the GOP) as Ax — 0. From R = kNZ, it is clear that this divergence is 
caused by the small level spacings s = {E\ — E^)/A which are distributed according to 
p{s) oc for small s [^]. Indeed / dss^~'^ converges near s = for (3 > 1 but diverges 
logarithmically for (3 = 1. The expression (|5lD for Z can be regularized by introducing a 
small spacing cut-off at e = 5. From Y2{e) ^ 1 — vr^e/G + 0{e''^) 0, we find 

Z, = -hn5 + 0{l) . (65) 

To regularize the perturbative expression (^) we examine the conditional average for | 
A7A P at fixed H{x) (see f^)) 



I A7A |2 |fixed//(x)= - 7 T \2 I I'' ' (66) 

^ ^^A ^f^) 

written in terms of the unfolded energies and the scaled parameter. This expression is well- 
defined for any given member H{x) of the GE ensemble at x. However, for a member with a 
level spacing that is too small, the approximation (0), which is of the order | Ax 1*^, breaks 
down unless we also make | Ax |^ smaller. Consequently, the average over all H{x) at fixed 
I Ax I'' leads to divergence. To avoid that, we use a small spacing cut-off S (as for Z above) 
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that depends on Ax. Noting that the parameter which determines the goodness of the 
approximation (|66D is | Ax I'' we make the conjecture that a regularized perturbative 
expansion can be obtained by assuming 5^ cx| Ax \^. This is equivalent to relating two 
(unfolded) energy scale S"^ oc (Ae)^ where the latter is associated with a change Ax in the 
scaled external parameter. This regularization leads to 

R = kN In I Ax I'' +C(1)) . (67) 

Our conjecture is supported by an exact result for N = 2 which exhibits this behavior (see 
Section VII), and is confirmed numerically for large N below. In analogy with the universal 
ratio R/D for the GUP (|63D we now have a universal ratio for the GOP (for small Ax) 

^ = In (Ax) + Const. (GOP) . (68) 



Eq. (pSD can be used to determine the scaling factor D (in cases it is immeasurable) from 
the measured R. 

Fig. 1^ (left panel) presents the scaling ratio R/D computed with a finite Ax as a 
function of In Ax. Shown are results of GP and Anderson model simulations compared to 
our conjecture for the orthogonal case (p^). Similar calculations are shown (right panel) for 
the unitary case as a function of Ax and they converge to (0) in the limit Ax — > 0. The 
agreement is very good. We remark that the smaller Ax is, the more simulations should 
be performed to get good statistics due to increasing fluctuations. Consequently, the GOP 
divergence would be difficult to observe experimentally. 



C. Perturbative Calculation of the Width Correlator 

We now calculate cr(x — x') for a symmetric dot perturbatively to leading order. A 
convenient way to calculate the width correlator is to use the relation 



cr(x-x') = 1- -i^^ , (69) 



where AF^ = Fa(x') — Fa(x). Using Fa(x') =| ■yx{x') p in ( ^3]) we get 



^^lll, + c.c] +... . (70) 



After averaging 



{r'x - rxr = E T^^^i^^^/^^^r, + c.c. (7i) 



over H{x') while keeping H(x) fixed, similarly to the calculation (^61) , we are left with 



2fl2 T 

(n - rxy = --(1 - f )E T^^^rr^r.F, (72) 

P f,^x (Ex - ^m) 
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where the average is over H{x). Eqs. (|T7| , |49| , 



531) then give 



(r^ - TxY = AT^TxR I Ax r ■ (73) 



Using the GE relation T^Tx = f ]^ (ri - Ta ) (valid for /i ^ A) together with (H) we 
obtain in the level velocity scaling 

R 

cr(x - x') = 1 - Pj^^ I X - x' I" . (74) 
For the GUP we use Eq. (|63D to find 

cr(x — x') ~ 1 — I X — x' 1^ (75) 

as —>■ oo. When we compare ( [75| ) with the leading-order behavior of the squared 
Lorentzian ( ^2]) for the GUP case we find a2 = v^/ti" ~ 0.55. The discrepancy with 



the least-square fit value quoted below (|3^) indicates that (p2|), while being a good approx- 
imation, is not exact. Indeed, including higher-order terms in the perturbative calculation 
introduces odd powers of | x — x' \'^^'^ into cr (see Section V.D). 

In the GOP the perturbative approximation ( |5B| ) diverges. Using a regularization similar 
to that invoked in the previous section amounts to replacing R in (|53|) by R as given by 
(|57|). Using ( |B5D this leads to the following short-distance behavior of the GOP conductance 
correlator 

cr(x — x') ~ 1 — — — - I X — x' I'' ( — — In I X — x' I +Const. j . (76) 



This non-analyticity near the origin indicates that the Lorentzian ( P^ which provides a very 
good approximation for Cr(x — x') is not the exact result. We remark that the N = 2 case 
exhibits the same behavior (see Section VII). 



D. The Rotation Scaling as the RMS Width Velocity 



In Ref. [|T3| we have shown that the quantity (Ae^)^, which measures the diffusion of 



the unfolded energy levels in x, is a universal function of (Ax)'' (even when the latter is not 
small). Similarly, one can consider the quantity (ArA)^/r;^ where AF^ = Fa(x') — Fa(x) 
as a measure of the diffusion of the "unfolded" widths Tx/Tx- Based on our discussion 
from Section UFA, this must also be a universal function of (Ax)''. We have computed this 
function for the GOP and GUP using the simple GP cosxiJi -|- smxH2, and it is shown in 
Fig. |. 

Using the perturbative expression ([73D, we find that for small Ax 



A^ ^ . . .R 



n^^^'^' = ^ivTW2^ ' '""^ ' ■ ^^^^ 

For the GUP we expect to find R/D R/D = 7r^/3 in the limit Ax 0, while for the 
GOP R/D should have the logarithmic universal behavior (|68D at small Ax. 
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Eq. ( ffTD leads us to an alternative definition of the rotation scaling 



1 1 




2 


ATI 


1 Ax 1 


n 



R - ^^rf^ • (78) 



In most physical applications, the parametric dependence is analytic so 77 = 2, and the scaling 
factor a/R is just the RMS of the rate of change of the normalized width as a function of 
the parameter 



We remark that the definition ([79| ) of the rotation scaling has the advantage that there 
is no phase ambiguity in T{x). Furthermore, in quantum dots experiments, the measurable 
quantity is the width (i.e. conductance for a symmetric dot), so ( [79| ) allows us to extract 
the scaled parameter directly from the data. 

It is possible to calculate perturbatively higher order terms in (|77D. In the GUP, while 
the leading order converges, the next order does not. For example, one of the terms of 
order Ax^'' includes a diverging GUE average ^/{E\ — E^Y. If we introduce a small 

level spacing cutoff 5, this average behaves as 1/5. Using the same regularization as in 
Section IV. B, i.e. 5 oc Ax^^^, we find a perturbative term of order Ax^''^^. Carrying our 
regularization procedure to higher order terms we make the following conjecture regarding 
the analytical behavior of (^) for finite Ax. For the GOP 



4fi 



2 

- =1 Ax 1^ [/ii(Ax) - /i2(Ax) In | Ax |^] , (80) 



where hi and h2 have power expansion in (Ax)''. For the GUP 



=1 Ax I'' [/ii(Ax) - /i2(Ax) I Ax I"/'] 



ATI 



(81) 



where hi and /i2 also have power expansion in (Ax)'' (but are different from the corresponding 
functions for the GOP). Our conjecture is supported by the exact results for the N = 2 GP 
derived in Section VII. Note that for small Ax the functions in brackets in (|80| ) and (^T]) are 
just R/D. 



E. Semiclassical Determination of the Rotation Scaling and the Correlation Field 

To obtain a universal behavior it is necessary (for rj = 2) to scale the external parameter 
by \/D or alternatively by a/R. These scaling factors are system-dependent and in principle 
could be determined semiclassically. In this section we present a nonrigorous yet simple 
semiclassical derivation of the rotation scaling factor R. The WKB wavefunction inside the 
dot is given by \1/ ~ Ae"^^^^, where S is the classical action. In the limit h ^ 0, the sensitivity 
of ^ to X is due to the x-dependence of the action and we find, from ( |5^ ) 
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In the important case where a; is a magnetic field B, it is possible to calculate 
explicitly. In particular we would like to derive an expression for the correlation field Be 
defined by Br = \fRB = B / Be- The change in the action when the field is changed by Ai? 
is given by 

AS =^ J AA -di. (83) 

In the following we shall distinguish between open and closed dots. For open dots we 
write AS" = ^-^AB, where G is the area (times 2%) swept by a classical trajectory. Using 
fB3) we then find 



R ~ 972$^ , (84) 

where $0 = hc/e is the fiux unit. In a ballistic chaotic open dot the area distribution is 
given by P(6) oc e~"'='l®l so that 6^ = 2/ali and R = l/a^^^g- The correlation field is 
then given by 

Be = aa^o . (85) 

We note that this is exactly the correlation field that appears in the semiclassical derivation 
of the conductance correlator in ballistic open dots We also note that oc t~^J'^ where 
Tesc is the mean escape time from the dot |3^ . 

We now turn to the case of closed dots (the Coulomb blockade regime) which is the 
one relevant to this paper. For such dots the electron does not escape but continues to 
accumulate phase until time r = h/A, where A is the mean level spacing Since 



F <^ A, the decay time is much longer than r, and it is that time scale r that determines 
the correlation field. The action difference AS" is then calculated for this time r. We write 

AS = - r dt ip{t) , (86) 
c Jo 



where (p{t) = AA(r(t)) ■ v(t). To calculate (A^*)^ we follow the method used in ||39| 

(ASy = 2— dt / dt' ip(t)ip(t + 1') ^ 2—T / dt <^(0)(^(t) . (87) 
Jo Jo & Jo 

Since the magnitude of the velocity v along the orbit is constant, the length i travelled in time 
t is given by £ = vt. Measuring as a function of I and for an electron with unit velocity 
(w = 1), we have dt ^(O^W = v dl (p{^)(p{i), where (p{i) = AA ■ v = \ABd{i) 
with d being the distance of the corresponding orbit segment from the origin. Defining a 
dimensionless geometrical constant ag [p9| 



ag = ^-3/2 / diO)d{£) , (88) 



we find from (R2D 



21 



R{^Bf = agn'rvA-^/^ (^^)' . (89) 
A^^'^/v the time for the balhstic electron at the Fermi energy to cross 



If we denote by 
the dot, then the correlation field can be written as pO 

1 /^£o 
Tc^/a^\ T A ' 



where 



r 



Tr,. 



V2mEA 



(90) 



(91) 



and we have used A = 2iTh'^ / niA. We note that our result (|90D is similar to the one obtained 
for the field that is required to break time-reversal symmetry [|^,^,^. From (^) and (|9l|) 
we find that BcA/ $o oc J\f~^/^ where J\f is the number of electrons in the dot. In the diffusive 
regime we expect the correlation field to be given by a similar expression except that the 
ballistic crossing time is replaced by the time it takes the electron to diffuse across the dot 
Tcr = ^h/Ec where Ec is the Thouless energy. We then find 



(92) 



in agreement with the fluctuation correlation field obtained in the supersymmetry method 
18| , and with the field strength that is required to break time reversal symmetry in disor- 




dered systems • 

Finally we remark that for open dots the semiclassical conductance correlator 



written in terms of the scaled field Br as c, 

form Cg = (l + A_B^/a2r) (see Eq. 
converted a2 of Section III.B to its value in the rotation scaling). 



^ can be 

g - (l + AB^y . For open dots we find a similar 
2f) for (3 = 2) but with ~ 1-15 (where we have 



V. UNIVERSAL DISTRIBUTIONS OF THE WIDTH VELOCITY 

In this section we restrict the discussion to the different iable parametric dependence 
(?7 = 2) where the factors associated with the level velocity scaling and the rotation scaling 
are related to the second moment of the level velocity and width velocity, respectively. One 
may be interested in higher moments of these quantities or in general in their distribution. 
In fact, for any finite Ax we can define an average velocity of the unfolded levels and widths 
with respect to the scaled parameter x 

Ax r Ax ^ ^ 

Their respective distributions Pax{v) and Pax(^) are symmetric and universal for any given 
Ax. We now apply the perturbation techniques of Section IV to derive expressions for these 
distributions in the limit of small Ax. 
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We first address the simpler case of the level velocity distribution P{v) = 
6{v - Ae/Ax) = (27r)-V^e*''*e-**^^/^^. In first order perturbation theory v = {H[^ - 
Hxx)/Ax, and the conditional average on H{x') at fixed H{x) can be easily done by notic- 
ing that H'y^^ — fHxx is a Gaussian variable with zero mean. Since cummulants higher than 
second vanish for a Gaussian variable, we find 



e" 



Ifixed Hi.)= e-^*^(^^/^^)^ = e-^*^ , (94) 



where we have used (Ae/Ax)^ = 1 + 0{Ax'^) at fixed H{x). Since the expression found in 
(0) is independent of H{x), the average over the latter does not change it and 

P{v) = (27r)-i/2g-it-2 _ (95) 

Thus in the limit of small Ax, the level velocity distribution is simply a Gaussian. In Fig. ^ 
(top panel) we show (for both the orthogonal and unitary cases) the universal distribution 
(P5|) (solid lines) and compare it with Pax{v) for Ax = 0.1 (circles), as calculated from the 
GP. Only at larger values Ax > 0.2 we observe deviations from the limiting distribution 
(p5|). In Fig. I (top panel) we show the universal form of the distributions P/^xiv) for such 
larger values of Ax. 

We now turn to the more interesting calculation of PAx{r) in the limit of small Ax. To 
leading order in perturbation theory (see ([70|)) r = (AxF)"^ J2 (^ e^-e 7a 7m + '^•'^•^ ^ 
Gaussian variable at fixed H(x) with zero conditional mean. Therefore 



e-*r- Ar/A. 1^^^^ ^-i.^f-(Ar/A.)^ ^ (gg) 

where the conditional average of the squared width's rate of change is 



V^^Vi r,F,/F^ 



Taking the Fourier transform of (p6D, we can express P(r) as a GE average over Gaussian 
distributions 



P(r) = - 
^ ' 2 



TT 



^^x (^M - ^a)2 



E 



-1/2 



4 V ^A'^A/r;^ 



(98) 



Thus we have reduced the process average in the original expression for P(r) to an ensemble 
average. Fig. ^ (bottom panel) shows this limiting distribution (|98D by a solid line, while 
the distribution Paz('") for Ax = 0.1 is shown by circles. Notice that the variance of the 
limiting distribution P(r) is just R/D = in the unitary case and R/D in the orthogonal 
case. In Fig. |^ (bottom panel) we show the universal distributions Pi^x{f) for larger values 
of Ax. We use a logarithmic scale in order to show more clearly the deviations ar larger 



values of r from the limiting form (|98|) . 

We observe from Fig. ^ that for both the level and width velocity distributions the 
deviation from the corresponding limiting form occurs faster (with increasing Ax) for the 
GUP than for the GOP. Also in each symmetry class (i.e. orthogonal or unitary) the change 
in the distribution as Ax increases is slower for the width velocity distribution than it is for 
the level velocity distribution. 
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VI. CONDUCTANCE CORRELATIONS IN ASYMMETRIC DOTS 



In present experiments in the Coulomb blockade regime the dots do not have reflection 
symmetry. In this section we calculate the conductance correlator and its analytic short- 
distance behavior for such asymmetric dots. We also express the rotation scaling in terms 
of RMS conductance velocity. 



A. Universal Correlators 

For dots without reflection symmetry F^, 7^ F^, and Eq. (y) should be used for the con- 
ductance G. If in addition the leads are asymmetric we need in general two different channel 
correlation matrices M^^,, MJ^, of dimensions to characterize the dot, assuming that 

channels at different leads are uncorrelated. However, Eq. ( p4D still holds for F''^ separately 
and the conductance correlator cg{x — x') depends only on the eigenvalues (w^)^''" of M''^^, as 
in the symmetric case, since the invariance of the GP under global unitary transformations 
implies that the choice of (f)^f does not affect averages involving {(f)^f \ ip). 

We first investigate cg{x — x') for an asymmetric dot with one-channel leads. This is 
the case which is relevant to present experiments in the Coulomb blockade regime where 
one channel dominates in each lead |2D[. The strongest deviation from the symmetric dot 



correlator is found for symmetric leads F' = F*". The correlator for this case is calculated 
from simulations of the GP (^Tj) and presented in the top panel of Fig. ^ (dashed) and is 
compared with the symmetric dot (|32|) (dotted). As in Section III.B we can also fit this 
correlator to (^) but now we find ai = 0.37 ± 0.04 for the case of conserved time reversal 
symmetry (GOP) and 0:2 = 0.54 ± 0.04 for broken time-reversal symmetry (GUP). Also 
shown in Fig. ^ are results from the Anderson model as described in Section III.B. In the 
top right panel (broken time- reversal symmetry), the parameter dependence is introduced by 
folding the cylinder into a torus and applying an additional magnetic field B' perpendicular 
to the torus plane. This amounts to setting 6ai3 = — $'/$o = ^ in (^) if fa, fp are neighbors 
along X and using 9 as the parameter. In the case of highly asymmetric leads, the small 
width can be neglected in the denominator of Eq. (^, and cq reduces again to the width 
correlator. We conclude that the conductance correlator for an asymmetric dot with single- 
channel leads varies between the two limiting cases of symmetric leads and a symmetric dot. 
We note that this variation is smaller in the unitary case. We also find that it takes a large 
asymmetry of the leads to see significant changes in the correlator. To study the effect of 
increased number of channels on the conductance correlator, we computed cg{x — x') for 
an asymmetric dot with multi-channels symmetric leads (M^ = M^) assuming uncorrelated 
and equivalent channels. The GP results for 8 channels leads and the Anderson model 
results with 4x4 leads are presented in Fig. ^ (bottom) and show cg{x — x') approaching 
the symmetric dot result. We emphasize that in this case the Anderson model results are 
not expected to agree with the GP since their channel correlation matrices differ, thus the 
agreement we find numerically suggests that the correlator in the asymmetric case converges 
to a limiting form as the number of channels increases and that this limit is the symmetric 
dot correlator. 
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B. Perturbative Expression for the Conductance Correlator 



Using the techniques of Section IV. C, we can calculate cg{x — x') perturbatively (to 
leading order in Ax) for an asymmetric dot with arbitrary leads. We use the notation 
G = ieyh){n/2kT)gx where gx = r^r^^/(r^ + Tl) and define Ag^ = g^ix') - gxix). For a 
small change in x we have 



I \ 2 



-pr 
^ X 



Ar 



A ' 



(99) 



where Tx = T'^^ + r\ is the total width of the level A. Since the two leads are uncorrelated 
we get 



r, 



1^ (Ar,)2 + (^Ia^ (Arl)2 



(100) 



r^^/r^ under an x-independent orthogonal (unitary) transformation, we find 

1 



Using the same methods of Section III.B, and taking advantage of the invariance of 



Lr 



r 

Ta 



Tr(M'''')2( ^ ) {ATliy 



(101) 



where r^'{ is the partial decay width to a single normalized channel in the corresponding 
lead and the last equality is valid in the limit of large A^. 

To calculate the expectation values in ( |101| ) we first perform the conditional average over 
H{x') (keeping H{x) fixed), for which we can use relation (|72[) to obtain 



AR 



Tr(M 



r\2 ^ X 



tI\2 / A 



r^,ri^ + Tr(M^ 

i A / VJ- A 



^ ix^ l^l 



(102) 



where /x denotes a level different from A. Using the relation cg{x—x') = \ — {AgxY/2{g\—g'^, 
and the independence of the eigenfunctions A and /i (in the limit oo), we obtain 



cg{x — x') 
where the constant 6/3 is given by 



R 



hp— I X - X' 



^' 1 2 



(103) 



9i 



Tr(M02(f^ + Tr(M')2(^ 



N{gl - 9l) 



(104) 



We remark that the quantities averaged over in ( |104|) are functions of the independent 
(normalized) partial widths f^'^ through F^'^ = N J2ci'^c^)^^cXj that each partial width 
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is distributed according to a distribution in (3 degrees of freedom with the same mean 
■^cA — 1/^- Therefore the constant bp is completely determined by the eigenvalues (Wc^)^ 
and the number of channels in each lead. 

In the important case of a dot with one channel in each lead and symmetric leads 
(f ' = f we find h = 7/4 and 62 = 3. It follows that for the GOP 

cg{x — x') ~ 1 — -{x — x'Y \ In I a; — x' I +Const\ , (105) 

4 V 6 ; 

while for the GUP 

CG(a;-x') ^ l-7r2(x-x')^ . (106) 



C. The Rotation Scaling as the RMS Conductance Velocity 

In asymmetric dots' experiments it is not possible to measure the width directly and 
it is therefore useful to express the rotation scaling directly in terms of the conductance 
velocity rather than the width velocity (see Section IV. D). This is accomplished through the 
perturbative expression (|102|) 



2Rbf,^^^ I Ax P 



9' 



(107) 



Eq. ( p.07| ) leads us to yet another definition of the rotation scaling 



1 /AG, ^2 



Gl \ Ax 



(108) 



where 



2b, 



In the case of one-channel symmetric leads, we find 



1 1 /AGy^^ 



^~ 7G\[ Ax 



(109) 



for the GOP, and 



5 1 dGx 
2AGI [ dx 



(110) 



for the GUP. 

In practice, one can calculate R from the measured conductance peaks according to 
( |109|) or ( |110|) , and then use the universal ratio R/D, as given by Eq. ( |68D or Eq. (|63D, 
to determine the diffusion scaling factor D for the GOP or GUP, respectively. The scaled 
parameter x that leads to universal correlations is now easily determined from (pOl). 
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VII. THE N = 2GP 



For = 2 it is possible to calculate analytically the width correlator discussed in this 
paper. While the results are different from the — > oo hmit, it is nevertheless instructive to 
study the N — 2 case. In particular the analytical properties of the rotation scahng factor 
and the width correlator are similar. 



A. Two-matrix distribution 

A GOE 2x2 matrix is real symmetric and thus has three independent real variables 
Hii,H22 and H12. We parametrize it in terms of its eigenvalues Ei,E2 and the rotation 
angle a that characterizes the orthogonal SO{2) matrix R that diagonahzes H. Introducing 
the level spacing E = Ei — E2 and average energy E — (Ei + i?2)/2, the GOE probability 
distribution is Vi{E, E, a) (x\ E \ er^^/"-'' e~^^/'^°^ , where -00 < E,E < 00 and -n/2 < a < 

Tl/2. 

To solve for various two-point correlators, it is sufficient to calculate the joint two-matrix 
distribution P{H, H') where H = H{x) and H' = H{x'). Introducing the quantities E, E, a 
and E', E' , a' for H and H' , respectively, we find 



P^{E,E',E,E',a,x) 



EE' 



4(7ra2)3(l - /2)3/2 



where x = ~ Q^' is the rotation angle of R ^R'. Introducing polar coordinates for E, E' 
through E = d^cosO, E' = ^sin^ and integrating over E, E' we obtain 

Pi(C,^,«,x) = 8(7ra2)4-/^) ^' ' """"^^ ' e-^^-^^''^^^-^^-)^^/^'^^^^-^^) , (112) 

where < C < 00, — tt < 6 < n and — 7r/2 < a,x < '^f'^- 

For wavcfTinction correlators that are independent of the spectra, we can integrate over 
( and 6. In order to keep the two levels at x and x' in the same order we require EE' > 0, 
and therefore the integration range for ^ is < ^ < 7r/2 and — tt < 6 < — 7r/2. We obtain 

2f1 - P) 

^i(0 = + 7^ \/t~D ' ^ = ^^«^2x . (113) 

A 2 X 2 GUE matrix is characterized by two real variables Hii,H22 and one complex 
variable H12. We parametrized it by its eigenvalues Ei,E2 and its diagonalizing unitary 
matrix 

/ e**^ cos ae*'^ e~*"^ sin «e*'^ \ (^^A\ 
" [ -e^'^sinae"*'^ e"*"^ cos ae"*'^ J ' ^ ^ 
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where O<a<7r,O<0< 27r, — vr < ip < it. We remark that the matrix H is independent of 
ip but it is more convenient to use the parametrization ( |1 14| ) since it describes an element of 
the group SU (2) whose invariant measure is D[U] =\ sin 2a \ dadipdil) . The GUE distribution 
is then given by V2{E, E,a,ip,(j)) oc E"^ sin2ae~'^^^ ^"-^ e~^^ ^"^"-^ . 

To parametrize the two-matrix distribution we introduce E, E, U and E', E', U' to de- 
scribe H and if', respectively. Instead of U and U' we then use U and U = U^^U'. Parametriz- 
ing U as in ( |114D but with angles ^ (replacing a, 0, ip) and using the invariance of the 
group measure D[U'] = D[U], we find 

P2iE,E',E,E',a,x)^ 

E^E'^ I sin2asin2x | e-2(^'+^''-2/-^^')V(i-/^)e-(^^+^''-2/i^^'™s2x)/2a2(i-/^) _ ^^^g) 

Integrating over E,E', 0,$ and ip,"^, and using polar coordinates (,9 in the E,E' plane, 
we find 

P2{C,0,a,x) oc C^sin22^ | sin2asin2x | e-^'(^-2/^''^2ecos2x)/2a2(i-/2) _ ^^^g^ 
If ^ and 6 are integrated out we have 

oc| sin 2a sin 2x | 92(0 
(72(0 = ^Y^^ f^e + ;^|^arctan^/|^j ; e = /cos2x. (117) 



B. The Width Correlator 

We take a channel *^ ~ ^ g ^ ^^^^ eigenfunction ip{x) = ipi{x) so that for the 

GOP r = Rf^ = cos^a and F' = Rf^ = cos^ a' = cos^(a; + x)- The correlation between F 
and F' is then given by 

1 — . . ^ 1 



FF' = cos"^ a cos^ x + ^^^^^ sin^ x " -^cos"^ a sin 2a sin 2% = -(1 + 2cos2 x) > (US) 



where we used cos^a = 3/8 and sin^ 2a = 1/2. Using F^ = 3/8, we obtain for the width 
correlator 



cr = cos2x. (119) 

For the GUP F =| Un p= cos^a and F'^ =| U[i p= cos^ a' =| cos a cos x — 
sin a sin xe"^**-''^^"^^ P, and therefore 



FF' = cos^ a cos^ x + ^^^^^ 2a sin^ x ~ -cos^ a sin 2a Ree M<P+'^) sin 2x 

= (l+^;^^)/6, (120) 



where we have used cos^a = 1/3 and sin^ 2a = 2/3. Since F^ = 1/3 we obtain again Eq. 
( |IT9|) for cr. 
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Eq. ( [119| ) can be calculated analytically for both classes. In the orthogonal case we have 
Cr =^ / cos2x\gi[0 - 9i{-i)\dx = — ; — / r. ^TTT^^X • (121) 



47r Jo ^ 4 7o (l-/2cos2 2x)=^/2 

The integral in ( |121| ) can be evaluated to give 



-/(l-/^)2Fi(3/2,3/2;2,/^), (122) 



where 2F1 is a hypergeometric function. In the unitary case we find 

_ (1 - p?" 



/ sin2xcos2x[5'2(0 - 5'2(-0]c?X 



(1-/2C0S2 2x)3/2 

Here too the integral can be done exactly to give 



13=2 o^/-'- ~ P I n'^f'^ " ^ ■ f /10/l^ 

— ' +2 — — arcsmj . (124) 



7lf 71 f 



C. The width Correlator at Small Distances 



We now examine the behavior of {AT)'^/AT'^ = (1 — cr)/2(/5 + 1) for small Ax = x — x' 
(see Eq. (^)) by expanding the exact solution for cr in 1 — / ^ Pn'^{AxY /8. 

For the GUP we find that ( |124| ) has a power expansion in (1 — fY^'^ and we obtain the 



general form ([81| ) with i] = 2 suggested in Section IV. D where 



^i(Ax) = ^ + ^(Ax)2 + ... 

MA.) = ^ + ^(A.)^ + .... (125) 

The leading order behavior agrees with the one derived from perturbation theory. We first 
use (|63| ) where for = 2 the relevant Z is Dyson's Zu = 1/4 (see Eq. (0)) to find 
R/D = From (^ we then obtain 



in agreement with the leading term in (|125|) . 

For the GOP we find after expanding ( |122| ) the general from (^) where 

, , 2 31n2-l-ln(7rV8) , 2, 

o4 

h2{x) = ^ + ^(Ax)^ + 0{Ax') . (127) 
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Thus the leading order terms in hi and /i2 are constants, which is consistent with the large 
perturbative result (|68|) for R/D (although with different constants). It is instructive to 
repeat the perturbative GOP calculation for the special case N = 2 and see whether our 
regularization procedure for this case yields the correct expansion. We obtain ([77|) with 
R/D given by (p2| ) but now 



4J e 8 V 16 y 



:i28) 



where we introduce a small level spacing cutoff S (in units of the mean spacing) and Ei is 



the exponential integral function. Using the identity Ei(— t) = C + \nt + J dt'{e 



where C is the Euler constant, we find 



C + ln 



7l'6 

16 



Making the assumption 5^ oc (Ax)^, the quantity (^) is then given by 

r2 



(ATT 

4r^ 



(Ax 



;^;^2 



TT 



^ In(Ax)^ + Const. 



(129) 



(130) 



in agreement with (|8lD and the leading order in (|127|) . 



VIII. CONCLUSION 

In this paper we discussed the parametric correlation function cg{x — x') of the con- 
ductance peaks in quantum dots in the regime of Coulomb blockade oscillations. Using 
the framework of the Gaussian process, we demonstrated its universality and obtained its 
functional form as well as its exact analytic behavior at short distances for both cases of 
conserved and broken time-reversal symmetry. For a symmetric dot we proved that cq is in- 
dependent of the channel configuration, while for asymmetric dots with multi-channel leads, 
our results indicate that the conductance correlator approaches its symmetric-dot form as 
the number of channels increases. We proposed a new scaling procedure, the rotation scal- 
ing, which uses a scaling factor that can be extracted directly from the conductance peaks 
and is therefore experimentally measurable. Our main result (^) can be readily tested 
experimentally. It would also be interesting to apply our rotation scaling to chaotic systems 
other than quantum dots. 

This work was supported in part by the Department of Energy Grant DE-FG02- 
91ER40608. We acknowledge CM. Marcus for useful discussions. 



APPENDIX A: USEFUL GE RELATIONS 

In this appendix we derive various useful partial width amplitudes correlations in the 
Gaussian ensembles. In the following 0i, 02 denote two orthogonal and normalized channels, 
while ijjx, ipij. ci-re two eigenfunction (A 7^ /i). jcX = {(pc I "^a) are the partial width amplitudes. 
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We first derive the distribution of tlie partial width amphtudes 7ia, 72A, 7i/i, 72/x, which 
is given by 



N N N N 

P(7ia,72A,7im,72m) = j n(^V^.A#,g5(E I V^^A I' -1)5(E I V^^M I' -l)5(EV'aV^^M) (Al) 

j=3 i=l i=l i=l 

Using two sets of spherical coordinates for ipj\ ; j = 3, . . . , N and ipj^ ; j = 3, . . . , N in 
^lTI) we find 

P(7)(x[(l-i?D(l-i?y-|5|f ^ , (A2) 

where Rj^ =\ 71 a P + | 72 a P and S' = 7x^71/^ + 72a72m- Front ([A2|) we can calculate various 
moments and correlations of the partial amplitudes. 

However, most of the useful relations can be derived in a simpler way. From the exact 
distribution of a single amplitude -P(7ia) oc (1— | 71A P)^^ ^, we find 

1 P + 2 



N 

By squaring the relation \ 'jix P= 1, taking an average and using (^) we find 



1=1 



1 13 



I 71A PI 72A P = j^j^ . (A4) 

N 

Similarly, if we average the square the relation J2 I 7ia P= 1 and use we obtain 

A=l 

hirn^=^^^- (A5) 

To derive the correlation between the widths of two eigenfunctions to decay into two orthog- 

7V 

onal channels, we use J2 I 7jA PI Ijfi P = 1 together with 



P{N-l) + 2 1 



I 7u PI 72. P = ^;v + 2 N{N-1) ■ 
Finally, by averaging the square of the modulus of the orthogonality relation of the two 

N 

eigenvectors J2 llxliii = 0, and using relation (|A5|), we find 

i=l 

1 1 ^ P 



7u7im72a72. = -^^1 7iA PI 7iM P = " ^^r + 2 iV(iV - 1) ' ^^^^ 
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APPENDIX B: PERTURBATIVE CALCULATION OF THE CROSS-CHANNEL 

AND EXCHANGE CORRELATORS 



To prove (|27|) to leading order in perturbation theory we consider the quantity 



{T[^ — Tlx) (^2x ~ ^2x) where the subscripts 1,2 correspond to an orthogonal pair of nor- 
malized channel vectors 0i^2, and Tix{x) =\ ■jixix) p=| {(pi \ i'xix)) p, i = 1,2. Using (|56D 
we have 

which upon averaging gives 

(r;,-ru)(r'2,-r2A) = 2 {Yixi2xii,i*2, + c.c) i? | Ax r . (B2) 

The cross-channel correlation is then calculated from ( |B2| ) using 



riA^A - r^A = riAr2A - r2A - (r^^^ - riA)(r2;^ - r2A)/2 . (B3) 

After also using the GE relation (o^?moV72^ + cc) /2 = f ^ (TW^ - T^) we get 



N 



VixT'^x - TiA r'2A = (riAr2A - TiA Tsa) ( 1 - Pj^R I Aa: I" ) . (B4) 



Comparing with the autochannel correlator ( JT^ ) and using the GE relation rir2— Fi r2 
j—^ {ViTi — Ei^^, we recover (|27| ) to order | Ax I''. 



N- 

Relation ( P^ ) can be similarly proven in perturbation theory. Writing the correlator of 

interest in the form (7u72a7u72A + c-c.) /2 = EiaEsa - | 7u72A - 71*a72a 1 72, we calculate 
the second term on the r.h.s. Here we have to distinguish between the GOP and the GUP. 
Eor the GOP 

2(7L7|: + 7ia72a7i/.72m)^I Aa;|\ (B5) 



7i*a72A - 7iA72A I' 



while for the GUP we find 



7u72A - iM2X ? = 4| 7iA 1^1 72;. ?R I Ax r . (B6) 



Using also the GOE relation EiaE2^ + 7ia72a7i/.72/. = ]v^-'^i^-'^2A and the GUE relation 
^l\^2^l = ^riAr2A, we obtain 

I {rixl2Xl[*xl'2x + cc.) = T^xTVx (l - f3j;^R I Ax . (B7) 

Comparing ( P7| ) with (0) and using riAr2A = f ]^ (j'lx " Tia^) , we recover pg| ) to leading 
order in I Ax 
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FIGURES 



FIG. 1. Universal form of the conductance correlator for a symmetric dot (|22|). Shown are GP 



simulation results (diamonds) obtained from (31) and their best fit to ( p2| ) (dashed). The circles 
are results from the Anderson model in a cylindrical geometry for leads of 2 x 2 point contacts, 
and with the parameter being the strength of an external potential. In the case of broken time 
reversal symmetry (right) a magnetic flux of ^/^o = 1/4 is applied. For comparison we also plot 
the GUP result on the left and the GOP on the right (dotted). The GP results are obtained using 
300 simulations with N = 150, taking the middle third of the states, while the Anderson model 
results are from the middle 200 eigenfunctions of the Hamiltonian (|33| ) with IF = 1 on a 27 x 27 
lattice. 



FIG. 2. GP (diamonds with statistical error) and Anderson model (circles) simulation results 
for the scaling factors ratio R/D computed with flnite Ax, compared to the prediction in the 
orthogonal ( |68| ) and unitary ( |63| ) cases. Notice that the ratio is plotted as a function of In Ax for 
the orthogonal case and as a function of Ax for the unitary case. Here x = \fDx is the diffusion 
scaling. The Anderson model results correspond to the cases shown in Fig. || but for 4x4 leads. 



FIG. 3. The universal form of (Ar)2/r function of {Axf for the GOP (left) and the 

GUP (right). The results are obtained from simulations of the process (^T]). The asymptotic value 
for large Ax is 4//3. 



FIG. 4. Top: the level velocity distribution P{v = de/dx) given by (^) (solid line). Shown 
by circles is the distribution P[v = Ae/ Ax) calculated for Ax = 0.1. Bottom: The width velocity 
distribution P(r = T~^dT/dx) given by (]98| ) (calculated for Ax = 0.01). Also shown is the 
distribution P{r = F^^AF/Ax) for Ax = 0.1 (circles). 



FIG. 5. Top: the average level velocity distribution P{v) (where v = Ae/Ax) for various 
values of Ax = 0.25 (squares), 0.5 (pluses), 0.75 (crosses) and 1 (diamonds). The solid line is 
the limiting Gaussian distribution (p5|). Bottom: the average width velocity distribution P(r) 
(where r = r~^Ar/Ax) for the same values of Ax as for P{v). The solid histogram is the limiting 
distribution (p8|). Notice that we use a logarithmic scale for P{r) in order to distinguish between 
the various cases at larger values of r. 



FIG. 6. Universal form of the conductance correlator (|Tj) in asymmetric dots. Top: GP 
simulation results with single-channel symmetric leads (dashed) and Anderson model results with 
1x1 leads (pluses and crosses). Bottom: GP simulation results with 8 equivalent channels in each 
lead (dashed) and Anderson model results with 4x4 point-contact leads (squares and circles). 
The Anderson model is as described in Fig. |l| where the parametric dependence is on the strength 
of an external potential, except in one case (crosses) where the parametric dependence is that 
of an additional external magnetic fleld in a toroidal geometry. For comparison we also plot the 
correlator for a symmetric dot (|3^) (dotted). 
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